f=inline('(x+1025/x)/2');
x0=3;
for i=1:20
    x0=f(x0);
    fprintf('i= %i,x0= %10.8f\n',i,x0);
end